Observation of the bottomonium ground state, rj b , at BaBar 
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We present the first observation of the bottomonium ground state J7(,(1S) in the photon energy spectrum using a 
sample of (109 ± 1) million of T(35) events recorded at the T(3S) energy with the BaBar detector at the PEP-II 
B factory at SLAC. A peak at E y = 921.21?, 'g(stat) ± 2.4(syst) MeV observed with a significance of 10 standard 
deviations in the photon energy spectrum is interpretated as being due to the radiative transition T(3S) — > 7»y(,(15). 
This photon energy corresponds to an r]i,(lS) mass of 9388.912 s (ate*) ± 2.7(syst) MeV/c 2 . The hyperfinc T(15)- 
rji,(lS) mass splitting is 71. 41^ (stat) ± 2.7(syst) MeV/c 2 . The branching fraction for this radiative T(3S) decay is 
obtained as (4.8 ± 0.5(stat) ± 1.2(syst)) X 10" 4 . 
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Bottomonium spectroscopy started thirty years ago with the discovery of the T(nS) resonances [1] . The spin-singlet 
i—, states hb(nP) and rjb{nS) have yet to be observed. In particular, the ground state of the bottomonium spectrum, 
776(15*), was still missing. The mass difference between the T(15) and the 7]b(lS), the hyperfinc splitting, is very 
important in understanding the role of spin-spin interaction in heavy quark bound systems and in testing calculations 
^"'i and predictions from various models such as Quark Models, pNRCQCD and Lattice QCD [2]. Predictions for the 
, hyperfine splitting vary from 36 to 100 MeV/c 2 [3]. 

We report on the observation of the bottomonium ground state r\b from the radiative transition T(3S) — > 777b [4]. 
1— I ■ Theoritical predictions for the branching fraction of the decay vary from 1 to 20xl0~ 4 [3]. The CLEO III experiment 
has published a 90% confidence level upper limit for the branching fraction B[T(3S) — > jrjb] < 4.3 x 10 -4 [5]. 
\ The data used in this study was recorded with the BaBar detector [6] at the PEP-II asymmetric-energy e + e~ 
storage rings. It consists of 28.0 fb 1 of integrated luminosity collected at a e + e CM energy of 10.355 GeV, 
corresponding to the mass of the T(35 l ) resonance. Samples of 2.4 fb -1 and 43.9 fb -1 recorded 30 MeV below 
the T(35) and 40 MeV below the T(45) resonances were used for background studies. The trajectories of charged 
particles are reconstructed using a combination of five layers of double-sided silicon strip detectors and a 40-layer 
drift chamber, all operated inside the 1.5-T magnetic field of a superconducting solenoid. Photons are detected using 
■ a CsI(Tl) electromagnetic calorimeter (EMC), which is also inside the coil. The energy resolution for photons varies 
from 2.9% (at 600 MeV) to 2.5% (at 1400 MeV). 
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2. BACKGROUNDS AND SIGNAL SELECTION 

The signal for T(35) — > jr/b is extracted from a binned maximum likelihood fit to the inclusive photon energy 
spectrum in the center of mass (CM) frame. The monochromatic photon from the decay will appear as a bump in 
the photon energy (E 7 ) distribution. For an rjb mass of 9.4 GeV, and the T(35) energy, the photon energy shall peak 
at 911 MeV. We are therefore looking for an enhancement in the E 1 distribution near 900 MeV. 

2.1. Background contributions to the E~ t distribution 

There are two main background contributions to the photon energy distribution The first contribution produces a 
smooth non-peaking background. It comes from continuum events (e + e~ — > qq where q = u,d, s, c) and bottomonium 
decays. The second contribution produces peaks in the E 1 spectrum, close to the expected signal position. It comes 
from two processes: 
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• The exclusive decay T(35) — > 7Xhj(2P); x&./(2P) — > 7T(15), J = 0,1,2. The second radiative transitions 
produce a broad peak centered at 760 MeV. As there are three transitions, we would expect to observe three 
peaks. However, due to the detector energy resolution and to the Doppler broadening, that arises from the 
motion of the Xbj(2P) states in the T(3S) CM frame, the three peaks merge into a single broad bump. 

• The radiative production of the T(1S) through initial state radiation (ISR): e+e~ — > jisrT(1S). This process 
produces a peak centered at 856 MeV. 

In order to extract the r\b signal, it is crucial to understand both the lincshapes and the yields of the two peaking 
background components. 

2.2. Signal Selection 

The selection criteria have been optimized by maximazing the figure of merit S/ y/~B, where S and B represent the 
expected yield for signal and background respectively. The signal sample is obtained from a detailed Monte Carlo 
(MC) simulation. There is no reliable event generator to model the various bottomonium decays. A small fraction 
(9%) of the data sample (in the region 0.85 < P 7 < 0.95 GeV) was used to model the background. In order to avoid 
any bias, this small data set was not used for the extraction of the signal in the final fit. 

As the i]b is expected to decay mainly through two gluons, one can expect a large track multiplicity in the final 
state. Events are selected by requiring at least four tracks in the event and that the ratio of the second to zeroth 
Fox- Wolfram moments [7] be less than 0.98. 

Photons arc first required to be isolated from all charged tracks, and their shapes are required to be consistent 
with an electromagnetic shower: the lateral moments [8] are required to be less than 0.55. In order to reduce the 
contribution from ISR events e + e~ — > jisrT(1S), candidate photons are required to be detected in the central region 
of the calorimeter —0.762 < cos(# 7i z,ab) < 0.890, where O^xab is the angle between the photon and the beam axis 
in the laboratory frame. 

We apply a cut on the angle 6t between the direction of the photon momemtum and the thrust axis [9] . The thrust 
axis is computed with all charged tracks and neutral calorimeter clusters in the event, excluding the photon candidate. 
Given that the r\b is a spin-zero resonance, the angle distribution for the signal should be flat. However, for continuum 
events the distribution should be peaking at the forward and backward directions. We require | cos#t| < 0.7. 

Finally we apply a veto to reduce photons coming from ir° decays. These photons represents the main source of 
background. A photon candidate combined with any photons in the event is required not to have an invariant mass 
within 15 MeV of the nominal ir° mass. The energy of the second photon in the tt° candidate is required to be larger 
than 50 MeV. 

These selection criteria lead to an efficiency of 37% and 6% on signal and background respectively. 
The optimization procedure was checked on data using the broad peak from the second radiative transition of the 
T(3S) — > 7%bj(2P);xbj(2P) — * jT(1S) process. It yielded to a very similar cut optimization. 

3. FITTING PROCEDURE 
3.1. Introduction 

The r]b signal is extracted, after all selection cuts are applied, using a binned maximum likelihood fit to the inclusive 
photon energy spectrum in the CM frame in the range 0.5 < E 1 < 1.1 GeV. 
There are four components to the fit: 

• non-peaking background; 

• Xf,j(2P) — > 7T(1S') peaking background; 
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Figure 1: The top figure shows the distribution from the T(4S) Off-Peak data sample. The bottom figure shows the 
background subtracted distribution. The 7isr¥(1S) peak is clearly visible. 

• JisrY (IS) peaking background; 

• t% signal. 



3.2. Probability density functions 

The non-peaking backround has been parametrized by the following probability density function: f(E 7 ) = 
A(C + exp[-aE 1 -/3E%]). 

As explained above, due to detector energy resolution and Doppler broadening, the three peaks from the x&j(2P) — > 
7 T(liS) transitions are merged. The three peaks have been modeled using a Gaussian modified with a power-law tail 
on the low side (Crystall Ball (CB) function [10]). The relatives rates and peak positions between the three peaks 
have been fixed from the PDG values [11]. For each Xbj(2P) lineshape, the parameters of the power-law tail have 
been fixed to a commom value. The PDF parameters have been determined from fitting non-peaking background 
subtracted the E 1 distribution where the signal region (840 to 960 MeV) has been excluded (see Figure 2). For the 




Figure 2: Fit to the background subtracted _E 7 distribution, with the signal excluded, for the determination of the Xbj{2P) — * 
yT(lS) PDF parameters. The blue dotted line shows the Xbj{2P) 7 """(IS 1 ) lineshape and the full purple line shows the 
7/SflT(lS) lineshape. 



final fit, all the Xb.j{2P) — > 7T(15) component parameters were fixed to the values obtained from this fit, except 
the yield. 

The PDF for the JisrT(1S) (i.e. ISR peak) peaking background component was modeled with a CB function. 
All the CB parameters were obtained from MC. Given the detector energy resolution, the ISR and signal peaks are 
likely to overlap. Depending on the mass of the rjb, the overlap could be large. In the final fit, it was therefore 
decided to fix the yield of the ISR peak. The rate of the ISR peak was estimated using data taken 40 MeV below 
the T(45) resonance (T(45) Off-Peak data). The top plot of Figure 1 shows the E 1 distribution for the T(4S) 
Off-Peak data after all cuts are applied. The bottom plot shows the same distribution after subtracting the non- 
peaking background. A clear ISR peak is seen. A fit with a CB functions yields to 35800 ± 1600 events. The yield is 
extrapolated to the T(3S) energy using the relative cross-sections, integrated luminosities and signal reconstruction 
efficiencies. The estimated jisrT(1S) yield is then 25200 ± 1700 events. The error includes systematic uncertainties. 
This is consistent with but more precise than the yield estimated with data taken below the Y(3S) resonance. 

The r\h PDF is modeled with a non-relativistic Breit-Wigner funtion (for the natural shape of the 775) convolved 
with a CB function which models the energy resolution. The CB paramaters were fixed from MC. MC experiments 
have shown that the width of the rjb had to be fixed in the final fit. It is not known, but theoritical prediction vary 
from 4 to 20 MeV [12]. We used a value of 10 MeV for the nominal fit. 

3.3. Fit to the full data sample 

The final fit to the E 1 distribution was performed with the PDFs described above. The free paramaters were the 
Xbj(2P) —> 7 """(IS) process yield, the non-peaking background parameters and the signal yield. Figure 3(a) shows 
the E 7 distribution and the fit result. In addition to the non-peaking background, only the Xbj(%P) — ¥ 7T(lS r ) 
broad peak is visible. Figure 3(b) shows the non-peaking background subtracted plot in the signal region. The 
Xbji^P) — ► 7 """(IS 1 ), jisrY(1S), and signal peaks are clearly visible. Figure 3(c) shows the background subtracted 
distribution overlaid with the fit result for the rjb PDF. 

The fitted rib signal yield is 19200± 2000±2100 events, where the first error is statistical and the second systematic. 
The systematic error has been obtained from varying the r\b width (from 5 to 20 MeV), the ^isnTiXS) yield within 
±1 a of the nominal value, and the PDF parameters within ±1 a. 

The rjb signal significance is estimated using the ratio log(L max /Lo), where L max and Lq are the likelihood values 
obtained from the nominal fit and from a fit with the rjb PDF removed, respectively. The significance the signal has 
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Figure 3: (a) Inclusive E 1 distribution with fit result. Only the Xbj(2P) — ► tT(1S) is visible, (b) Non-peaking background 
subtracted plot, with PDFs for Xftj(2P) — > 7T(15) peak (solid), ISR peak (dot), r)b signal (dash) and the sum of all three 
(solid), (c) All-background subtracted distribution. 



been conservativaly estimated with the following method: a fit to the data has been performed with all parameters 
entering the systematic errors moved by 1 standard deviation in the direction of smallest significance. This method 
yields to a 10 standard deviations significance. 



4. RESULTS 

The fitted r]b signal position is 917.41^8 MeV. A photon energy calibration shift of 3.8 ± 2.0 MeV is then applied, 
ft is obtained by comparing the fitted position of the Xbj(2P) — * 7T(15) peak to the PDG value. Applying the 
energy calibration shift, we obtain for the peak position of the r\b signal: E~ f = 921.2+ 2 ' 8 ^ 2-4 MeV. 

This yields to the % mass: M(r ?6 ) = 9388.91^3 ± 2.7 McV/c 2 . Using the PDG value of 9460.3 ± 0.3 MeV/c 2 for 
the T(15) mass, we determine the T(lS)-rjb mass splitting to be 71.4± 2 .-3 ± 2 .7 MeV/c 2 . 

Using the signal reconstruction efficieny and the number of T(3S) events, we estimate the T(35) — > 7 775 branching 
fraction to be (4.8 ± 0.5 ± 1.2) x 10 -4 , where the first uncertainty is statistical and the second systematic. The 
main systematic uncertainty is from the efficiency. We have compared the reconstruction efficiency of the \bj (2P) — > 
7T(15) peak between data and MC, giving a 18% error. 



5. CONCLUSION 

In conlusion, we have made the first observation of the bottomonium ground state, the r\b- The new state has been 
observed in the radiative decay of the T(35). The rjb is the most likely interpretation of the signal, although other 
hypothesis are not excluded. The mass of the rjb is 9388. 9+I3 i 2-7 MeV/c 2 , which corresponds to a mass splitting 
between the T(IS') and the rjb of 71.4+ 3 ' 1 ±2.7 McV/c 2 . The estimated branching fraction of the decay T(35) — > 7776 
is found to be (4.8 ± 0.5 ± 1.2) x 1CT 4 . 
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